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SUMMARY 


Calculations are made to determine the trajectories of liquid drop- 
lets introduced into the air disturbances generated by an airplane engaged 
in aerial spraying. The effects of such factors as the position at which 
droplets are ejected into the disturbances, airplane lift coefficient, 
and altitude are investigated. The distribution of deposit on the ground 
is computed for several droplet- size spectra, variations in the rate at 
which mass is ejected along the span, and lateral flight-path spacings. 
Consideration is then given to the problem of adjusting these factors 
with the aim of improving the uniformity and increasing the effective 
width of the deposit. 

The results indicate that the lateral dispersion of droplets is 
increased when the spanwise position at which particles are ejected is 
moved toward the wing tip. Greater dispersion also results when the air- 
plane lift coefficient or altitude is increased. 

With the spray discharged continuously along the span at a constant 
rate, the deposit has a maximum concentration at the plane of symmetry, 
and rapidly diminishing concentration beyond the wing tip. In such cases, 
it was found that the \miformlty and effective width of the swath co\ild 
be in 5 )roved by increasing the mass efflux rate with spanwise distance 
from the center line of the wing. When this is done, the lateral dis- 
tance between adjacent flight paths required for a given degree of uni- 
formity shows an increase amounting to from 0.4 to 1.0 semispan. Of the 
two droplet-size spectra considered, the one having the smaller mean 
diameter (200 microns) gives a more uniform deposit for a given flight- 
path spacing and the degree of uniformity is less sensitive to changes 
in the spacing of adjacent passes of the airplane. 


INTRODUCTION 


Althou^ airplanes have been used to disperse liquid sprays and 
dust for insect control and agriciilt'ural purposes for many years, only 
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recently has an appreciable amount of study been devoted to the problem 
of improving the efficiency and effectiveness of the operation. One of 
the . problems recognized as being iii 5 )ortant is that of obtaining a suita- 
ble distribution of the deposit. Experimental studies made by the Forest 
Insect Laboratory of Beltsville, Maryland (ref. l) have indicated that, 
in general, aerial spray deposits exhibit excessive peaks of concentra- 
tion beneath the airplane with diminishing intensity toward the outer edge 
of the swath. Although it is possible to improve these measured deposit 
characteristics by a tedious trial-and-error process, it is believed that 
a theoretical study of the in^jortant physical factors involved would pro- 
vide a better mderstanding of the mechanism which determines the lateral 
dispersion of particles ejected into the flow field behind an airplane. 
Such a study may also suggest possible methods of Improving the deposit 
characteristics which might otherwise be overlooked. 

In the present paper the paths of liquid spray droplets issued into 
the flow field behind an airplane are calculated. Considered in these 
computations are the effects of such factors as the droplet-size spectrum, 
spanwise location of nozzles, height of airplane above the groxmd, and 
aerodynamic flow field in the wake. Distributions of deposit are deter- 
mined for several variations of these parameters, with a view to inproving 
the uniformity and the effective width of the swath. The airplane used 
as an illustration in the analysis is the agrictiltxiral prototype desig- 
nated the Ag-1. 


SYMBOLS 


A wing aspect ratio 

b wing semispan, ft 

C concentration of deposit, Ib/Unit area 

Cp drag coefficient of droplet 

Cp lift coefficient of airplane 

D drag force on droplet, lb 

5 diameter of droplet, microns (l micron = 3*28 X 10“^ ft) 

Bjj mean diameter of droplets, microns 

g acceleration due to gravity, ft/sec^ 

m mass of droplet, sliigs 
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P = 

Q 

1 

R 

% 

r 

t 

U 

V 






x,y,z 


- y, 


0 

mass ejected per unit length along flight path, lb /Unit length 
velocity induced by isolated rectilinear vortex, ft/sec 
Reynolds nxunber of droplet 

Reynolds nvutiber of droplet moving through still air at 
velocity U 

radial distance from vortex center, ft 
time, sec 

flight velocity, ft/sec 
percent of total volume 

velocity con^jonents of air in y- and z-direction, 
respectively, ft/sec 

velocity con5)onents of droplet in y- and z-direction, 
respectively, ft /sec 

longitudinal, lateral, and vertical coordinates of droplet, 
respectively, semispans 

vertical coordinate of droplet with origin shifted to the 
bovind vortex, semispans 


oCjjRt 

a dimensionless parameter, 

24Ru 

r circulation, ft^/sec 

e dispersion of droplet diameters, microns 

lateral and vertical coordinates of trailing vortex cores, 
semispans 

^ji absolute viscosity of air, slugs/ft-sec 

V kinematic viscosity of air, ft^/sec 
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density of air, slvigs/cu ft 
density of droplets, slugs /cu ft 

o dimensionless parameter, 18 — — 

5 P^ 

T dimensionless time, Ut/b 

Subscripts : 

0 initial value 

g intersection of droplet path with ground plane 

y,z lateral and vertical directions, respectively 

t terminal velocity of droplet 

c correction due to bound vortex 


ANALYSIS 


The paths taken by spray droplets issued from an airplane are influ- 
enced by the flow field induced in the region of air throvigh which the 
airplane has flown. When drag forces acting on the droplets are pre- 
dominant, as in the case of small particles, the trajectories follow the 
streamlines of flow more closely than do the trajectories of larger par- 
ticles for which inertia forces are predominant. In aerial sprays, the 
smallest particles may remain in suspension in the air indefinitely, 
whereas the largest particles fall in paths practically uninfluenced by 
the surroxmding air flow. In this section the method of determining tra- 
jectories of aerial-spray droplets and the restating distribution of 
deposit is presented. 

Calculations of a similar nattire have been performed by various 
investigators of aircraft icing problems. In these studies the trajec- 
tories of small water droplets moving in an airstream past various bodies 
such as airfoils and cylinders have been calculated to determine area and 
rate of droplet in^jingement on the surface. Glauert, in reference 2, 
considered cases wherein the air velocities and droplet diameters were 
sufficiently small (R < l) so that the drag coiad be assumed proportional 
to the velocity of the droplet relative to the air (Stokes ' law of resist- 
ance). Later, Langmuir and Blodgett (ref. 3) extended this work and con- 
sidered the variation of drag beyond the limits of velocity and diameter 
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in which Stokes ' law is valid. It was fovind, after a consideration of 
the range of droplet Reynolds nxraibers encompassed in aerial spraying, 
that the deviation of drag from Stokes ' law should be included in the 
present calcxilations . 


Equations of Motion 

The motion of spray droplets in a space-fixed plane located behind 
the airplane perpendicular to the line of flight will be considered. A 
rectangular-coordinate system is used where vertical z distances axe 
measured in airplane semispans from the ground plane (positive upward) 
and lateral y distances in semispans from the line of flight. If v^ 

and w^ are, respectively, the lateral and vertical components of drop- 
let velocity and Vg^ and w^^ the corresponding coniponents of air veloc- 
ity in the plane, the equations of motion which establish dynamic equi- 
librium of the forces acting on a droplet are 


m 


dvd 

dt 


= % 


( 1 ) 


dwj 

m = Dz - mg (2) 

dt 


The resultant aerodynamic drag force on the droplet is 


D = 





2 


and the lateral con^ionent of drag which appears in equation (l) may be 
expressed 
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According to Stokes' law (see ref. k, p. 598) the drag of spheres at very 
low Reynolds nxmibers (R < l) is 


%tokes 






2 


Now, since 


R 


^ Pa^U^a - Vd)^ + (va - ^d)' 


the ratio of the actual drag to the drag based on Stokes' law becomes 

D _ *^dP 
^Stokes 24 


Since 


m 



% - 


PaSU 

P 


and 


T = 


b 


the nondimens ional equation of motion in the y-direction, after appro- 
priate substitutions into equation (l), becomes 


d 



dr 


_o. 

Ru W~\\j " uy 


( 5 ) 
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where 


a 


18 


b Pa 

^ Pd 


In a similar manner the equation of motion 
becomes 

<^(^d/U ) ^ _g_ CdR^ _ 

dT Ry 

Except for the incliosion of the term due to gravity, bg/u2, which 
appears in equation (4), equations ( 5 ) and (4) are of the same form as 
the equations of motion given in reference 5 . Since the left-hand sides 
of equations ( 5 ) and (4) represent the accelerations of the droplet in 
the lateral and vertical directions, a double integration of these equa- 
tions with prescribed initial values yg, Zq, (^dA^)Q> (wd/u)o 

defines the s\ibsequent positions of the droplet. 

In carrying out the integration of equations (l) and (2) it is 
necessary to evaluate R and Ci)R/24 as the particle progresses along 

its trajectory. The Reynolds nxmiber of the droplet is computed from the 
relation 



and from reference 5 ibe term CyR/24 is tabulated as a function of R 
in table I. 


in the z-direction (eq. (2)) 


Uj -y2 


(4) 


Air Velocities 

The air velocities v^ ajid Wg^ as functions of droplet position 

and time are needed to solve equations ( 5 ) to ( 5 ). The two most signifi- 
cant distvirbances behind an airplane are contributed by the wing and the 
propeller. Of these, the effect of propeller wash is largely localized 
to the vicinity behind the fuselage and its air distvirbance velocities 
are approximately in the direction of the flight path. Becaxise the 
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inclusion of these and other air velocities in the direction of flight 
would greatly complicate the analysis, and since such disturbances are 
believed to have a small effect on the lateral displacement of droplets, 
they have been neglected. The wing, on the other hand, sheds a sheet of 
vortices which soon roll up into two discrete vortex cores that are 
initially spaced somewhat less than a wing span apart. These so-called 
trailing vortices and the velocities induced by them in the yz-plane 
persist for a considerable time after the airplane's passage nniS appear 
to be of fundamental in^jortance in determining the lateral distribution 
of droplets which impinge on the groiind. 


Idealized flow pattern .- An idealized flow pattern is used in the 
present analysis to approximate the velocities which are induced behind 
a wing. The trailing vortices and their images reflected from the ground 
plane are represented by a system of four parallel rectilinear vortices 
which pass perpendicularly through the yz-plane at points A, B, B', and A' 
(see fig. 1) and extend to infinity ahead of and behind this plane. The 
vortices are of equal strength and have the directions of rotation indi- 
cated by the figure. The coordinates of the vortex cores of this system 
change with time, since the coordinates of one of the vortex cores moves 
at the velocity induced at that point by the remaining three vortices. 

As a resxilt, the trailing vortices are displaced downward and outward as 
depicted by the dashed lines in figure 1. 


If the equations which express the velocities of the cores and the 
coordinates of their paths (ref. 4, pp. 225 to 224) are combined and 
integrated with respect to time, the positions of the cores as functions 
of time are obtained. If a rectangular spanwise loading is assumed, the 
strength of the trailing vortices is expressed as 


CrUb 
r = -i— 
A 


and these relations then become 
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where t) ftriti ^ are, respectively, lateral and vertical coordinates 
of the vortex core in semispans and the consteoits of integration are 



The velocity components v^/U and induced by the vortex 

system at a point having coordinates y and z semispans are, by the 
Biot-Savart theorem. 


va C 


U 2rtA 


^ + z 




(t, + z)2 + (t) - y)2 (5 + z)2 + (ti + y)' 


t - z 


^ - z 


- z)2 + (tj + y)2 - z)2 + (ti - y)‘ 


(8) 


^a _ Cl 
U 2rtA 


T - y 


Ti + y 


(C + z)^ + (ti - y)^ + z)^ + (ti + y)^ 


T] + y T] - y 

(5 - z)2 + (tj + y)2 _ z)2 + (tj - y)^ 


( 9 ) 
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The streamlines indicating the direction of air velocity vectors at a 
given time are shown in figure 2. The effect of foliage and other 
obstructions to air flow near the ground is, of course, neglected by 
the idealized flow pattern. 

Beca\xse the induced velocities given by equations (8) and ( 9 ) become 
infinite at the position of the line vortices, it is necessary to con- 
sider the TnlnlTmim radial distance from these points at which these equa- 
tions remain applicable. In reference 5 approximate expression is 
derived for the radius of the rolled-up vortex cores behind a wing. It 
is shown by simple energy considerations that the core radius for a wing 
having an elliptical spanwise loading is 0.155 semispani In order to 
avoid the region of the vortex core, the spanwise starting positions of 
the droplets should be at least this distance from the center of the core. 

Velocities in a real fluid . - Since the determination of velocities 
in the assumed flow field is based on an ideal fluid, it is of interest 
to consider briefly the effects of viscosity on the flow velocities in a 
real fluid, or, more specifically, the dissipation of vorticity with time. 

To illustrate the effect of the decay of vorticity on the velocities 
induced in a real fluid, consider an isolated rectilinear vortex with an 
initial strength of Fq. At a subsequent time t, the circulation in a 

circle of radius r about the vortex (see ref. 4, p. 592 ) is 

r = Fo(l - (10) 


Hence the tangential velocity on the circle becomes 

Viscous = (n) 


and therefore 


^viscous _ K ^-r2/4vt 
'lideal ^ 

Equation (l2) is plotted against r for various values of time in figxire 3* 
The figure shows the extremely slow rate at which induced velocities a 
few feet from the vortex line are diminished as a result of viscosity. 
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An additional viscous effect not accounted for by the idealized flow 
pattern is evidenced in the form of a boundary layer of reduced air 
velocities adjacent to the ground. 


Trajectory Computations 

The sample airplane for which the trajectory computations were made 
is the Ag- 1 , a single -engine, low-wing, agric\iltiu:al airplane of moderate 
size. The liquid droplets are assumed to have the density of kerosene, a 
frequently used solvent for insecticides. The physical characteristics 
pertinent to these computations are given in table II. 

A step-by-step integration of equations (l) and ( 2 ) was accomplished on 
the Bell Telephone Laboratories X- 66 'jhh relay con^uter at the Langley 
Laboratory. It was found convenient to tabulate the variation of CpR ^ 2 k 

with R for several values of R and use a linear interpolation between 
the tabiilated values. The coordinates of the trailing vortex cores were 
determined in a similar manner by tabulating tj and ^ at suitably 
chosen values of t. The air velocities were con^juted at each point 
along the trajectory by evaluating equations (8) and (9) directly. 

A series of trajectories was computed for several flight conditions 
with selected values of the parameters a and Ry and with assigned 
initial conditions y^, Zq, and ^Wd/UjQ. In order to reduce the 

number of computed points along a trajectory and yet retain a reasonable 
degree of accuracy, the time increments diiring an integration were varied) 
the criterion used to select these time intervals was the change in the drop- 
let acceleration evaluated at two adjacent points. The integration was 
stopped when the particle intersected the ground plane or when its tra- 
jectory looped about the vortex core. 

The vortex cores are assumed initially to be at the height of the 
airplane and separated laterally a distance of 1 span. Particles are 
given starting positions at various points along a line connecting the 
trailing vortices, and the initial droplet velocity is taken to be the 
steady-state velocity of the droplet if it were falling under the influ- 
ence of gravity in a uniform airstream having velocity components equal 
to the induced velocities of the idealized flow pattern at the point of 
release. This assumption was made, rather than the assumption that the 
particles were released from rest at the trailing edge, in order to com- 
pensate somewhat for the error caused by neglecting the bound vortex in 
regions close to the wing (see appendix A). Actually, in the neighbor- 
hood of the wing, air velocities are induced by a bound vortex at the 
wing and, in addition, a continuous sheet of vorticity extending down- 
stream from the bound vortex. In the present case, where a sin5)lified. 
uniform spanwise loading is considered, the flow in this region may be 
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represented by a single horseshoe-vortex arrangement with a bound vortex 
at the quarter-chord point connected with trailing vortices of equal . 
strength which terminate at the wing tip and extend downstream to infin- 
ity. By neglecting the bound vortex, the initial downwash is underesti- 
mated, althoiigh this error in downwash is to some extent reduced by the 
forward-extending vortices assianed in the idealized flow pattern. 

At a distance of 2.0 semispans behind the bound vortex the velocities 
induced by the idealized flow pattern were fovind to be within 6 percent of 
the velocities induced by the horseshoe -vortex system. The approximate 
error in the vertical position of a droplet 2.0 semispans downstream of 
the wing, caused by an underestimation of the downwash ahead of this 
point together with the cou 5 )ensating assuirption of a finite downward 
ejection velocity, is evaluated in appendix A. 

A summary of the results of droplet-trajectory computations is given 
in table III. Here, the quantities yQ, yg, . v^g, w^^, and tg are 

given for the various droplet sizes and airplane flight conditions con- 
sidered in the analysis. The calculated droplet trajectories for these 
cases are plotted in figure 4. 


Distribution of Deposit 


From the foregoing trajectory conputation it is possible to determine 
the distribution of deposit for various assumed nozzle characteristics and 
spanwise rates of efflux. Consider first the case where a spray, conposed 
of particles of various sizes, is introduced into the flow field at point 
^yQ,ZQ^ and at the rate of Q pounds per unit length along the flight 


path. For each droplet size there is a corresponding position of impinge- 
ment on the ground therefore, droplets of size 5 and 6 + fall 


at yg and yg + respectively. The rate at which mass is deposited 

in element ZVg is Q Av where Av is the percentage of, Q which is 

represented by droplets having diameters between 5 and 6 + /55. The 
concentration of deposit at position yg is then 


C = Q 


Av 

Ay, 


g 


(15) 


which, in the limit, can be written 


c = 

<35 dy. 


(14) 
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Equation (lU) applies for a single nozzle, and by superposition the 
concentration of deposit can be determined for any number of nozzles. 

For a large ntunber of nozzles closely spaced along the wing, however, it 
is more convenient to consider mass as being ejected continuously along 
the span rather than at a finite number of discrete points. 

The rate of change of concentration at yg with respect to Q is, 
by equation (l^), 

^ ^ dS 

dQ dB dyg 

but, since Q is a continuous function of yQ, 


dQ 


dQ 

dyo 


dyo 


Therefore, the total concentration at yg due to a continuous mass flow 
along the span becomes 


C = 


0 ^ ^ ^ 
dyo dS dyg ^0 


(16) 


In order to con^jute the concentration of deposit by means of equa- 
tion (1^) or (16), suitable expressions for the terms dv/dB and dB/dyg 
must be determined. The relation dv/dB is defined by the nozzle or 
spray characteristics, and dS/dyg is associated with the trajectories 

of droplets. 


Nozzle and Spray Characteristics 

Reference 6 shows that the breakup of liquid jets, ejected into a 
high-velocity alrstream, is dependent on the Reynolds number of the mean 
droplet size and the relative viscosities of the liquid particles and the 
surrounding air. The nozzle bore was foimd to have no effect on the mean 
droplet size except at very low ejection velocities, in which case the 
droplets were more xmiform and had diameters approximately twice that 
of the nozzle. 
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An empirical expression for the mean diameter, which is applicable 
for all conditions likely to be encovmtered in aerial spraying, is given 
in reference 6. In the notation of the present paper the expression is 


5m = 


500vct 


1/5 


(17) 


where the mean diameter 5jjj is the droplet size in cm about which the 

largest percentage of the total volume is concentrated, V is the rela- 
tive velocity in cm/sec between the liquid jet and the air, and is 

the kinematic viscosity of the liquid in cn^/sec. 

The spectra of droplet sizes measured in reference 6 were symmetrical, 
so the droplet size at 50 percent volume was taken to be Moreover, 

the measiu-ed distribution of droplet sizes appeared to have the form of 
the so-called normal distribution found in the theory of probability. 

The equation of the normal distribution (see ref. 7) then provides the 
following mathematical e3q>ression for the term dv/dB appearing in equa- 
tions (ik) and (16): 


^ = 0-6T3 

dB e^2jt 


-1/2 


0-675(5 - Sm) 


(18) 


The term e in equation (18) is known as the "probable error" and 
in the present problem represents the dispersion or spread of droplet 
sizes about Bjjjj that is, e is such that 50 percent of the volume is 

composed of droplets having diameters between 5^-6 and Bjjj + e. The 

magnitude of e is then a measure of the viniformity of droplet sizes in 
a liquid spray having a mean diameter B^. 

Two sanqole spectra of droplet sizes are considered in the present 
analysis; these have mean diameters of 200 and 500 microns and a disper- 
sion e of 50 microns. The percentage of the total mass in terms of 
droplet diameter may be determined by integrating equation (18), a plot 
of which is given in figure 5 for the two droplet spectra under considera- 
tion. (See tabulated values of the probability integral in ref, 8. ) 
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Mathematical Representation of Trajectory Data 

An en5)irical relation which expresses the droplet diameter as a func- 
tion of yq and y-g for a given flight condition is derived in appen- 
dix B by applying the method of least squares to the trajectory data pre- 
sented in table III. An equation which was found to give a satisfactory 
approximation of the con^juted data is 


6 = + kgP + 


(19) 


where 



and the coefficients k^^, ^2, and k^ are functions of yq. It is seen 

that equation (19) satisfies the condition that the droplets which descend 
in practically vertical paths ^y^ Y(^ are very large, whereas the drop- 

lets which are carried away in the vortex flow ^y^_ ^ 00^ have finite diam- 
eters equal to or smaller than k^* 

Figure It indicates that as the droplet diameter is reduced the influ- 
ence of disturbance velocities becomes more pronounced, until finally a 
critical condition is reached wherein the path of the droplet encircles the 
vortex core. It appears from the limited nvmiber of trajectory computations 
available that once the particle encircles a vortex core it continues to do 
so for a considerable time, during which it is carried laterally outward 
with the trailing -vortex core. Eventually the induced velocities become 
dissipated sufficiently so that they can no longer sustain the particle 
and it falls to the ground; however, the mass represented by such parti- 
cles and also the mass "lost" throvigh evaporation is neglected in deter- 
mining the concentration of deposit. 

To determine the percentage of mass carried away by the vortex veloc- 
ities it is necessary to consider the droplet spectrum and the spanwise 
position from which droplets were ejected. Since coefficient k^^ in 

equation (19) represents the minimum droplet diameter which is recovered, 
in terms of the starting position, the percentage of mass neglected by 
the analysis is simply the percent volume corresponding to 6 = k^^ in 

figure 5 ' In figure 6 , the ratio of the mass recovered to the total mass 
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ejected is plotted as a fvmction of the droplet starting position for the 
two san5)le droplet-size spectra with = 1.2 and Zq = 0.5 semispan. 

When adjacent swaths are overlapped in order to cover a given area more 
uniformly, the part of the mass which is neglected by the analysis, but 
does not evaporate, is probably distributed more or less uniformly over 
the entire area. 


RESULTS AND DISCUSSION 

Distribution of Deposit With Discrete Nozzles 


When liquid is ejected behind the wing from a discrete point, such 
as a nozzle, the res\ilting concentration of deposit on the ground may be 
determined by evaluating equation (l4) with the aid of equations (18) 
and (19) qt various positions on the ground. The distribution of 

deposit so computed is presented in figures 7 to 9 for a unit mass efflux 
rate. These figures are intended to illustrate the effect of spanwise 
location of the nozzle, mean diameter of the droplets, airplane lift 
coefficient, and airplane height. 

The effect of the mean diameter of the droplets on the shape of 
the deposit ciirves is shown in figure 7 for the case of Cl = 1.2 and 

Zq = 0.5 semispan. Here, the deposit curves are characterized by high 

peaks of concentration when the nozzles are located near the plane of 
symmetry ^yq = 0.25 semispanj, whereas for nozzle locations in the vicin- 
ity of the wing tip ^y^ = 0.75 semispan) the concentration of deposit is 

less intense and the distribution is more uniform. This effect would be 
expected since the magnitude and curvature of the air-flow velocities are 
increased as the wing tip is approached, and, as a consequence, droplets 
of different diameters are dispersed more in this region than at points 
farther inboard. 

Curves for various values of mean diameter of the droplets indicate 
that as 5„, is reduced the deposit is distributed over a larger area, 
the center of which is displaced slightly farther outboard. 

The effect of airplane lift coefficient on the deposit curves is 
illustrated in figure 8, where the altitude is again taken to be 0.5 semi- 
span and 6m = 50O microns. At Cl = 2.2, the lift coefficient corre- 
sponding to the flaps-down condition, the deposit curve is of similar 
shape but displaced farther outboard than for the flaps -up case (Cl = I.2). 

It is of interest to note that the velocity of in5)ingement for the flaps- 
down case (see table TIl) is generally sli^tly higher than for the corre- 
sponding flaps-up case. The higher droplet velocities and also the stronger 
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air-disturbance velocities which accompany the \ise of flaps may provide 
greater penetration of the crop foliage by the sprayed liquid. 

A conqparison of the deposit curves calculated for flight altitudes 
of 0.5 and 1.0 semispan is shown in figure 9^ where Cl = 1*2 and 

6jjj = 300 microns. The figure indicates that the uniformity and lateral 
displacement of the deposit is increased with airplane altitude. 


Distribution of Deposit With Continuous Spanwise Mass Flow 

In practice, aerial sprays are generally issued from many nozzles 
along the span. In order to calculate the total concentration of deposit 
when there is a multitude of closely spaced nozzles, it is convenient to 
use equation (16), which is the integral form of equation (l4). Variations 
of dQ/dyo were assvuned and the resvilting concentrations of deposit were 

determined by integrating equation (16) nvonerically with the use of 
Simpson's rule. In performing these con5)utations the available trajec- 
tory data (table III) were extrapolated by means of equation (19) so as 
to encompass a range of spanwise starting positions from yo = 0 to 
Yq = 0.8 semispans. In the results which follow. Cl = 1.2, zq = O .5 
semispan, and the mass effl\ix rate is assmed to be zero beyond yo = 0.8 
semispan. 

In figure 10 are presented deposit ciirves evaluated for a uniform 
spanwise effl\xx rate (solid lines in the figure) and also deposit curves 
in which dQ/dyo was modified in an attempt to improve the uniformity of 

deposit (dashed lines). For the case of a uniform spanwise mass efflux 
rate the deposit curves for both of the droplet-size spectra considered 
show maximum concentration at the plane of symmetry. The deposit curve 
for 5jji = 500 microns (fig. 10(b)) is the less uniform of the two; it 

has a higher concentration at the plane of symmetry but at yg = 2.0 semi- 

spans it has about one-half the concentration of the deposit curve for 
6jjj = 200 microns (fig- 10(c)). The deposit is, of course, reduced to 

zero in all cases as yg becomes large. 

For the purpose of determining a modified spanwise efflux rate which 
wo\ild improve the \iniformity of deposit, assume that all of the droplets 
issued from the wing were of equal size (8 = Constant) . Each of the drop- 
lets ejected at a point yQ would then fall at a single point -y on the 

O 

ground in accordance with the trajectory characteristics of the droplet 
size in question. The concentration of deposit in this case would be 

C = 

dyo dyg 


( 20 ) 
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Therefore, the variation of dQ/dyo with Yq required to provide a 
uniform deposit of a constant-droplet-size spray becomes 


dQ _ Constant 


( 21 ) 


where dyQ/dyg is expressed as a function of Yq alone by the simulta- 
neous solution of equation (l9)^ with 6 equal to a constant, and the 
relation 


dyp _ 

<3yg 8&/8yo 


Although this procedure for determining the optimum spanwise variation 
of mass flow rate is true only for sprays of constant droplet size, the 
xmiformity of deposits having a variable droplet size shoixld be iniproved 
by varying the mass flow rate in accordance with eqviation (2l), taking 
the droplet size to be the mean diameter of the droplets. 

The modified effl\nc rate I in figure 10(a) is equation (2l) evaluated 
with 6 = 500 microns. The constant in equation (2l) was selected to give 
the same total mass flow per vinit length along the flight path as was assumei 
for the case of a constant spanwise flow rate (Q = 0.8 pound per \mit length 
The resvilting deposit is shown by the appropriate coded curve in figure 10(b 
Here, the peak deposit which occurs at yg = 1.1 semispans indicates that 

too large a portion of the mass is transported laterally outward. One 
method of overcoming this difficulty would be to evaluate equation (21 ) 
for a droplet diameter somewhat larger than 6 jq. 

The mass flow rate evaluated to give a uniform deposit with a constant- 
diameter spray of 350 microns is denoted by II in figure 10(a)., The corre- 
sponding deposit for the case of a spectrm of 300-mlcron mean diameter is 
shown in figure 10(b). Here, the distribution of deposit is seen to be 
approximately viniform up to yg = 1.2 semispans, and beyond this point it 

breaks off abruptly as do the other deposit curves in figure 10 (b); how- 
ever, the lateral extent of the deposits of variable mass flow rate is 
somewhat greater than the deposit calculated for a viniform mass flow rate. 

Since a more favorable deposit results when the diameter used to 
evaluate equation (21 ) is somewhat larger than the mean diameter of the 
droplet spectrvun, efflux rate I was used to modify the deposit curve for 
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5jjj = 200 microns. The resulting deposit is shown in figure 10 (c). The 
concentration of the modified deposit for = 200 microns is less 
intense but is uniform over a larger region than the modified deposit 
with 6jn = 500 microns. Moreover, there is no abrupt reduction in the 
intensity of concentration beyond yg = 1.2 semispans as was the case 

for 6m = 500 microns. 

In addition to varying the mass flow rate along the wing, there 
remains the possibility of regulating the mean diameter of the droplet 
spectrum at various nozzle locations along the wing. As indicated in 
equation (17), this covild be acconiplished throiigh control of the rela- 
tive jet velocity, which is dependent on the direction of the jet rela- 
tive to the surrounding airstream and the pressure under which liquid is 
ejected. 

investigate one such case, it was assumed that the mean diameter 
increased linearly from 200 microns at the plane of symmetry to 500 microns 
at yo = 0.8 semispan and that the spanwise efflux rate was held constant. 

In this example, the spectrum of smallest mean diameter (6 = 200 microns) 
was issued from the plane of symmetry in order to reduce the peak deposit 
which occurred in this region when sprays of similar spectra were issuqd 
at a uniform rate across the span. The distribution of deposit for this 
case is plotted in figure 11. It is seen that this deposit curve has 
characteristics similar to the modified deposit of figure 10 (b), where the 
iiniformity of deposit was improved by increasing the mass flow rate with 
distance from the midspan. 


Effect of Flight-Path Spacing on Distribution of Deposit 

Thus far, an attempt has been made to improve the uniformity of the 
lateral deposit of an isolated swath. In actual practice, however, uni- 
formity of the deposit is dependent upon not only the distribution of 
concentration within a single swath, but also the manner in which the 
swaths are superin5>osed, that is, the lateral spacing of adjacent passes 
of the airplane. 

Since fli^t-path spacing is an important factor with regard to 
economy of operation, its effect on the total distribution of deposit 
has been determined for several of the individual swaths shown in fig- 
ure 10. These results are presented in figure 12. In order to determine 
the total concentration of deposit, the ordinates of overlapping swaths 
are added as illustrated in figure 12(a). If an optimum deposit is taken 
to be one in which the minimum required concentration is distributed 
uni 'Pnr' Tnl y on the ground, the concentration which exceeds this required 
minimum value will constitute an excess of deposit. A ratio of the excess 
deposit (hatched area in fig. 12(a)) to the total deposit may therefore 
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be considered a measure of the ■uniformity of the svath. The variation 
of this ratio with lateral flight-path spacing is plotted in figure 12(b) 
for 6jjj = 500 microns and in figure 12(c) for 5jjj = 200 microns. The 

coded lines correspond to the efflux rates given by figure 10 (a). 

.These plots indicate that for a similar degree of xmiformity in the 
swath, the flight-path spacing can be increased by modifying the efflux 
rate in the manner sho'wn in figure 10(a). For example, if 10 percent 
excess deposit is considered tolerable, the modified flow rate would 
permit an increase in the fli^t-path spacing of approximately 0.4 semi- 
span for 6m = 300 microns and 1.0 semispans for 6m = 200 microns. A 

second conclusion to be drawn from figure 12 is that, of the two spray 
spectra considered, the spray with 6m = 200 microns provides a more 
vtniform swath for equivalent flight-path spacings. Furthermore, the 
uniformity of this swath appears to be somewhat less sensitive to flight- 
path spacing. This consideration would be of importance in cases where 
it is difficiilt to estimate the optimum spacing between adjacent passes 
of the airplane. 


CONCLUSIONS 


Resvilts of the foregoing analysis of factors affecting the dispersion 
of spray droplets issued from a moderate-size agricultural airplane indi- 
cate that: 

1. Greater lateral dispersion of aerial sprays will result when the 
position at which droplets are ejected is moved toward the •wing tip, when 
the mean diameter of the droplet spectrum is reduced, or when either the 
airplane lift coefficient or the altitude is increased. 

2. The uniformity and effective width of the swath is iniproved by 
increasing the mass effl\ix rate with distance from the midspan. 

3 . As con^ared with a spray spectrum of 300 -micron mean diameter, 
the spray spectrum with a mean diameter of 200 microns allows greater 
flight-path spacing and produces a more uniform deposit. Moreover, the 
degree of uniformity is less sensitive to changes in the spacing of 
adjacent passes of the airplane. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 4, 1955* 
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APPENDIX A 


EVALUATION OF EE®OR IN DROPLET POSITION DUE TO PROXBUTY 

OF THE BOUND VORTEX 


When a particle is released from the trailing edge of a wing, the 
velocity induced by the boxind vortex may represent a substantial part of 
the total velocity acting on the particle during the initial part of its 
motion. In the present analysis an attempt to account for this addi- 
tional downwash was made by simply assuming that particles were ejected 
downward from the trailing edge. The initial velocities of the particles 
were assumed to be equal to their steady velocities when falling throvigh 
a uniform stream having an air velocity equivalent to the velocity induced 
by the idealized flow pattern at the Initial position of the droplet. In 
this appendix, the errors introduced by such an assumption are evaluated 
for the case in which the wing is represented by a horseshoe vortex. 

In computing these errors it is convenient to determine an alr- 
velocity correction which, when added to the idealized flow velocities, 
will provide the velocities induced by a horseshoe-vortex system. The 
required correction is seen to be the velocity Induced by a horseshoe 
vortex of span 2 b with legs extending forward to infinity in the direc- 
tion of flight. The strength of this horseshoe vortex is the same as 
the strength of the vortices in the idealized flow field arid its sense 
of rotation is such that, when superimposed upon the idealized vortex 
arrangement, the forward -extending legs cancel the vortices which would 
otherwise extend to infinity ahead of the airplane. The correction veloc- 
ity induced by the horseshoe vortex reflected from the ground plane is 
small in relation to the velocity induced by the bound vortex at the wing 
and in the present approximation will be neglected. 

For convenience the origin of coordinates in this case is taken at 
the center of the bound vortex with x positive in a rearward direction 
arw^ z' positive in an upward direction. The downwash-velocity correction 
(denoted by the subscript c) in the xy-plane is given by the relation . 
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At X = 2.0 semispans, the velocity on the mldspan given by equation (A1) 
is approximately 5.8 percent of the downwash induced by the idealized 
flow field when t = 0 and the airplane altitude is 0.5 semispan. Con- 
sequently, the error in droplet position associated with the error in 
induced downwash beyond x = 2.0 semispans will be assumed to be neg- 
ligible. Moreover, it will also be assmed in the following calculations 
that the vortex cores are held fixed during the time required for the air- 
plane to travel 2.0 semispans, and that the downwash- velocity correction 
does not change with the departure of the droplet from the xy-plane. The 
last two assumptions lead to overestimation of the magnitude of the down- 
wash correction whereas the first leads to underestimation of its duration. 

If the equation for vertical motion of a droplet (eq. (2)) in the 
idealized flow field is subtracted from a similar equation in which the 
downwash-velocity correction (eq. (Al) ) and the increment in drop veloc- 
ity due to this downwash correction have been included, the resulting 
equation when is constant is 


dr 


§>1 

=1 

II 

^ 0 


Ry 2k 

0 

c; 

1 


(A2) 


vhicli may be written 


d^z*g ^ dz'g _ 

Ru2J^ dr Ru24\U/c 


(A3) 


An integration of equation (A3) over a distance extending from the 
trailing edge of the wing (x = 0.3 semispan) to a position 2.0 semispans 
downstream of the bound vortex gives the approximate vertical droplet- 
position error caused by using the idealized flow pattern rather than a 
horseshoe vortex in the neighborhood of the wing. If the horizontal 
droplet position (in semispans) downstream of the bound vortex is given 

by Jit + 0 . 3 , where U is the airplane flight speed, the additional 

downwash acting on a droplet released from the trailing edge 

may be determined by substituting t = x - 0.3 into equation (Al), where 

T = Hi,. It was found that over the interval 0.3"^ x ^ 2.0 semispans 
b 

and for y values equal to or less than O .75 semispan, equation (Al) 
was essentially independent of spanwise distance y. An en 5 )irical relation 
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which gives a satisfactory approximation of equation (A1) over this 
interval is 


kitk 

Cl 



c 


= -4.09 + 6.60t - 2.62t^ 


(A4) 


where 0 t ^ 1.7 

The drag term CpR/24 of a droplet falling at a steady-state con- 
dition in a uniform stream is the same as if the droplet were falling at 
its terminal velocity thro\igh still air. Therefore, equation (4) for a 
steady condition reduces to 


Cp^t ^ M 
24 o Xj2 


(A3) 


where the subscript t denotes a terminal velocity condition and this 
velocity is in a negative direction. Equation (A3) may be written in a 
more convenient form as 


cpRt ^ ^ 

24 ° u2 


(a6) 


which, together with the functional relation between CpR/24 and R 
(table l), defines CpR-t/24, the constant coefficient in equation (A3). 


Equation (A3) may now be put into the form 


d2z'c dz'c 

dr2 <iT 



+ a^^T + 



(A7) 


where 
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The general solution of equation (A7) is 


z'c = Cie-^''" + (ao 




(A8) 


With the initial conditions 


T ^ 0 


c 

z' = £. = 0 

dT 


the arbitrary constants become 


C, - -C = ^ ^ ^ 

Cl---C2-a _2^,5 


(A9) 


a a-' 


When eqtiation (a 8), with constants defined by equation (A9) > is 
evalxiated at t = 1-7, the value obtained is the approximate error in 
vertical droplet position caused by using the idealized flow rather than 
the flow induced by a horseshoe vortex. The downwash is \anderestdjnated 
when the idealized flow field is used, so the correction is additive in 
a downward direction. 

It remains to determine the influence of the ass\m 5 )tion \ised to 
con 5 )ensate approximately for neglect of the bound vortex, namely, that 
droplets are ejected at a finite downward velocity 


. % (mo) 

.uLo y \Vt 


Let Az' be the total difference in distance traveled by a particle 
ejected into an airstream at a finite velocity ([eq. (AIO)) and a similar 
particle released from rest into the same airstream. If the airstream is 
considered to have a uniform velocity Wa,/U and the drag term Ci)Rt/24 
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is again taken to be constant, (in semispans) from the differential 

eq\jation of droplet motion becomes 


/ 



(All) 


The approximate total error in vertical droplet position which may 
be expected as a result of the assvuned starting conditions used in the 
trajectory computations is then the difference between equations (a8) 

Ani^ (All) . This error was eveiluated at several spanwise starting posi- 
tions for the flight conditions = 1.2 and Zq = O .5 semispan and 

for droplet diameters of 200 and 5 OO microns (the mean diameters of spray 
spectra to be considered) • 

The calculated errors in vertical droplet position (in semispans) 
for these conditions are as follows: 


yo^ 

semispans 

Error* (semispans) for - 

6 = 200 microns 

6 = 500 microns 

0.25 

-0.017 

o.ooi<- 

.50 

-.011 

.015 

• 75 

■ 

0 

0 

.055 


*A negative sign indicates that the error is in 
an upward direction. 


Similar computations were made for the error in lateral droplet 
position caused by neglect of the horseshoe vortex. The corrected 
lateral droplet position was inboard of the calculated position; however, 
the magnitude of the correction was found to be insignificant. 


/ 




26 


NACA TN 5032 




APPENDIX B 

LEAST-SQUARES REPRESENTATION OF TRAJECTORY DATA 


In the con5)utation of the distribution of deposit for the case of a 
discrete nozzle the unknown coefficients ki, kp, and k^ and the span- 

wise starting position yq in equation (19) are constant. It is desired 
to determine the values of kq, kp, and kj in the formula 

6 = kj^ + kgP + kjP^ (b1) 


which give the best approximation of n data points (table III) at the 
partic\ilar yo value in question. The principle of least squares pro- 
vides a method of con^juting these coefficients which minimizes the sum 
Of the squares of the differences between the values of 6 in table III 
and the values of 8 con5>uted from equation (B1) . Thus the three unknown 
coefficients may be determined by solving the following set of simultaneous 
equations : 


Y^8 = kqn + k25IP + ^5X^P^ 

5Ip8 = kiI~P + k2YIp2 + k5YIp5 > 

YL = kq H p2 + k2 H P5 + k3 P**- 


(B2) 


The deposit for discrete nozzles was evaluated only at the yQ posi- 
tions for which trajectory computations were available. This limitation 
obviated the necessity of interpolating to determine the coefficients at 
other yo positions; however, in determining the deposit for the case of 
a continuous spanwise distribution of mass flow (eq. (16)), the variations 
of kq, 1^, and k^ with Yq must be known. From experience gained 

while determining the coefficients of equation (B2) for various values of 
Yq, it was fo\md that the variation of these coefficients with yq could 

be satisfactorily approximated by the relations 


= ^2X0^ + ^3X0^ 
k2 = b^Yo + b2Yo^ I 

^3 = = 1^0 + ^ 270 ^ 


(B3) 
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Equation (B1) then becomes 

6 = a2yo^ + ^3^0^ + ^lyo^ + i32yo^p + °iyo^ <^2yo^^ 


Equation (b 4) contains six unknown coefficients, and for the case in 
which Cl = 1.2 and zq = O.5 semispan (the flight condition for which 
the deposit with a continuous mass flow was evaluated) there is a total 
of 20 data points (see table III). Applying the method of least squares 
to equation (Bk) provides a set of six simultaneous equations: 



= a2^ + 

yo^ + ^i5Z ^ 0 ^ + 

t>2 Y. ^ 0 ^ ^ 


Cl 5: p2yo5 

+ C2 5 :; P%o^ 



= a2 ^ yo5 + 


^2Y ^ 0 ^ + 


^ P^yo^ 

+ C2 ^ 



= 82 Pyo^ 

^^Y ^ 0 ^ bi'^p2yo2 

+ h^Y p^yo^ + 


^lYl ■* 

- C 2 II p5yo^ 





+ h2Y ^^yo^ + 


ci^ p5yo5 H 

■ C 2 XI p^yo"" 



a2 p2yo5 + 

a^Y ^^0^ + ^1^ ^^ 0 ' 

- + b2^ P^O^ ■*■ 


ciUp^yo^ + 

C2y~ P^yn^ 
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^ ^0^ = 32^ ^3 Y. ^0^ ^2^ P^yO^ + 

ciY" i^yo^ + ^2 51 p^yo^ 


These equations may be solved for the xmknown coefficients (a2^ 
and so forth) by any convenient method of solving linear simultaneoixs 
algebraic equations. 

The following numerical values of the coefficients in equation (b 4) 
were determined for the case in which Cj, = 1.2 and Zq = O .5 semispan: 


= 

575.7 

aj = 

- 157 . 1 ^ 

hi = 

216.9 

b2 = 

-49.76 

Cl = 

-0.0557 

C 2 = 

-8.l4l 


The con?)arison of the points determined by trajectory computations eind 
the enpirical representation of these points by equation (b 1<-) is shown 
in figure 13 . 
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TABLE I.- VALUES OE CdR/24 AS A FUNCTION OF R* 



TABLE II.- AIRPLANE AND SPRAY CHARACTERISTICS 


Wing semispan, b, ft 19. 5 

Wing aspect ratio, A 5.55 

Density of droplets, p^, slugs/cu ft I .55 

Density of air, pg^, slugs/cu ft 0.002575 

Lift coefficient of airplane. Cl 1.2; 2.0 

Flight velocity, U, ft/sec 85.0;62.6 
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TABLE III.- TABULATED RESULTS OF DROPLE 7 T-TRAJECTORY DATA 
(a) Cl = 1-2; Zq = 0-5 semispan 


semispans 

yg. 

semispans 

(*) 

5 , 

microns 

ft/sec 

Wdg. 

ft/sec 

sec 

0.25 

0.506 

210 

2.09 

-1-97 

2.88 

•25 

Ml 

275 

2.08 

-3.08 

2.26 

•25 

.585 

375 

2.06 

-4.21 

1.71 

.25 

.347 ■ 

500 

2.01 

-5.65 

1-31 

.25 

.315 

700 

1.85 

-8.75 

- 

• 95 

.50 

1-675 

150 

10.58 

-1.11 

5.05 

.50 

1.130 

210 

9-95 

-2.12 

2.18 

.50 

• 925 

275 

7.84 

-3.03 

1-77 

.50 

• 795 

375 

6.45 

-4.40 

1.42 

.50 

.710 

500 

5.15 

-6.83 

1.15 

.50 

.628 

700 

4.66 

- -9.16 

.84 

.625 

(**) 

150 



— 

.625 

1.670 

210 

8.76 

-1.67 

2.42 

.625 

1.250 

275 

11.47 

-2.91 

1.61 

.625 

1.006' 

375 

10.05 

-4.44 

1-35 

.625 

• 895 

500 

8.93 

-5-71- 

1.00 

.625 

• 787 

700 

6.77 

-9.51 

• 75 

• 75 

(**) 

210 



— 

• 75 

1.77 

275 

6.56 

-2.56 

2.26 

.75 

1-35 

375 

11.05 

-5.84 

1-25 

• 75 

1.10 

500 

11.01 

-6.02 

• 95 

• 75 

• 95 

700 

8.82 

-9.66 

.65 


The y„ values are without cross wind; To account for a 

o 


cross wind, the product of tg and the cross-wind velocity com- 
ponent may be applied as a correction to the tabulated yg values 
**Con5>utations were discontinued before droplet intersected 
ground plane. ^naca]^. 
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TABLE If I.- TABULATED RESULTS OF DROPLET-TRAJECTORY DATA - Continued 

(b) Cl = 2 . 2 ; Zq = 0.5 semlspan 


yo^ 

semispans 

yg^ 

semispans 

(*) 

5 , 

microns 

Vdg, 

ft/sec 

^dg. 

ft/sec 

sec 

0.25 

0.692 

150 

2.56 

-1.36 

3-,72 

.25 

• 556 

210 

2.71 

-2.06 

2.61 

.25 

.h8o 

275 

2.76 

-2.98 

2.03 

/ .25 

' -^5 

375 

2.76 

-4.26 

1.54 

.25 

• 570 

500 

2.72 

\- 5-97 

1.23 


.350 

700 

2.46 

-8.91 

^.94 

.50 

(**) 

150 


/• 



.50 

1.290 

210 

15.14 

-2.45 

1.86 

.50 

1.025 

275 

11.51 

-3.16 

1.52 

.50 

• 855 , 

375 

8.79 

-4.53 

1.25 

.50 

.750 

500 

7.64 

-5.67 

1.02 

.50 

.653 

700 

6.26 

- 9.64 

•76 

.625 

2.305: 

210 

6.16 

- 1*93 

5*51 

.625 

1.11-30 

275 

15.00 

- 2.68 

1.45 

. .625 

1.095 

575 

15.90 

-4.54 

1.06 

.625 

. 9^0 

500 

11.60 

-6.57 

•87 

.625 

.815 

700 

9.11 

- 10.10 

.67 

•75 

(^) 

. 210 



~ ^ - 

.75 

2.210 

, 275 

. 5-94 

- 2.42 

3.11 

.75 

1-^75 

575 

14.03 

- 3-59 

1.18 

.75 

1.155 

500 

14.92 

-6.23 

.75 

•75 

.985 

700 

12.64 

- 10.40 

.65 


* 

The jTg values are without cross wind. To account for a 
cross wind, the product of tg and the cross-wind velocity com- 
ponent may be applied as a correction to the tabxilated yg values. 




Computations were discontinued before droplet intersected 
ground plane. • 
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TABLE III.- TABULATED RESULTS OF DROPLET-TRAJECTORY DATA - Concluded 

(c) Cl = 1.2 ; zq = 1.0 semispan 


yo^ 

semi spans 

yg, 

semispans 

(*) 

5, 

microns 

ft /sec 

^dg’ 

ft /sec 

tg, 

sec 

0.25 

1.750 

150 

6.56 

-1.13 

6.55 

.25 

.851 

210 

5.56 

-2.12 

4.47 

.25 

.500 

575 

3.01 

-4.21 

2.55 

.25 

.360 

700 

1.97 

-8.81 

1.75 

.50 

(**) 

150 

— 


— 

.50 

1.925 

210 

i«-.l6 

-1.90 

5.60 

.50 

1.070 

375 

6.29 

-4.08 

, 2.74 

.50 

.7^0 

700 

4.21 

-8.58 

1.70 


*The y-g values are without cross wind. To acco\mt for a 
cross wind, the product of tg and the cross-wind velocity com- 
ponent may he applied as a correction to the tabulated yg values. 

**Con5>utations were discontinued before droplet intersected 
ground plane . 
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Figure 7*- Effect of the mean droplet diameter on distribution of deposit 
for discrete nozzles. Cr = 1.2: zn = 0-5 semispan. 



Figure 8.- Effect of airplane lift coefficient on distribution of deposit 
for discrete nozzles, zq =0.5 semispan = 500 microns . 





Figure 9 -- Effect of airplsine height on distribution of deposit for discrete 

nozzles. Ct = 1.2; &jn = 500 microns. 
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Sh/ath no. 2 
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(a) Overlapping of adjacent swaths. 
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(c) bm = 200 microns. 


Illustration of the excess deposit of overlapping swaths and 
its variation with flight-path spacing, 




